Hydrologic Design 4501

Chapter 3: Solar Radiation and Surface Energy Balance

Ardeshir Ebtehaj
University of Minnesota
Table of Contents

1- Radiative Energy

Heat energy can be transferred via different mechanisms including:
  1. Conduction: Transfer of heat energy (or any conserved quantity like mass) through molecular diffusion. This is the main mechanism in transfer of heat energy in solid materials.
  2. Advection: Transfer of a heat energy via bulk displacement of matter, which can take place in fluids or soft solids.
  3. Convection: Transfer of a heat via bulk displacement of matter and turbulent diffusion. Convection, can take place in fluids and soft solids or mixtures where particles can move past each other due to turbulence. Turbulent diffusion is much stronger than molecular diffusion and is due to random and chaotic swirls of motions of fluid particles. When a flow is turbulent, molecular diffusion is negligible compared to the turbulent diffusion. In turbulent flows, the fluid velocity is often high and thus the trasport mechnisam is due to a mixture of advection and turbulent diffusion, which we often call convection.
  4. Radiation:Transfer of heat energy by means of electromagnetic waves, with or without any intervening physical medium.
The Sun is the main source of heat energy behind hydrologic, atmospheric, and oceanic circulation. Because of the differential radiative heating and thus temperature gradient from equators to poles, heat fluxes are formed to balance the uneven distribution of heat.
As we discussed, this differential heating generates large-scale pressure gradients and thus circulation patterns from equators to pole and vice versa that the drive the hydrologic cycle. The circulation cells (e.g., Hadly or Ferrel cells) are often called convective cells because convection is the main mechanism behind the large-scale transport of heat and water vapor in the Earth's atmosphere.
Electromagnetic waves are formed by elementary particle, called photons, travel at the speed of light at different wavelengths λ [m]:
where is the speed of light, T is the wave period . Often, we express , where is the frequency in Hertz [Hz]. As shown, the wavelength is the distance between two crests of a wave.
A body of mass with a temperature above zero Kelvin emits radiative energy, which will be characterize later through the Planck's law.
The electromagnetic energy (EM) spreads across different wavelengths forming a spectrum shown in the following figure.
The wavelength of interests for hydrologists may include ultra-violet, visible, infrared, and microwave bands.
It is worth nothing that almost 50% of the Earth heating is coming from the Sun's radiative energy in infrared bands, which has important implications on global warming.
To properly quantify the EM radiative energy we first need to define some definitions and terminology as follows:
Radiance: The radiant energy coming from a specific direction passing through a unit area perpendicular to that direction per unit time , where steradian is the unit of solid angle capturing the directinality of EM waves. Spectral radiance is the radiance per wavelength or frequency .
Irradiance: The radiant energy passing through a unit area per unit time integrated over all solid angles . Similarly, spectral irradiance is the irradiance per unit wavelength or frequency
Schematic of the solid angle and its unit (steradian). A steradian is the solid angle facing an area of . Clearly, a sphere is steradians.

1-2 Blackbody and Graybody Radiation

A blackbody is a hypothetical physical object that absorbs all incident EM radiation.
Planck's law: The emitted radiation by a blackbody at temperature T [K], in terms of spectral radiances, as a function of wavelength λ or frequency ν is expressed as follows:
where is the Boltzmann constant and represents the Planck's constant.
Planck's law as a function of temperature and wavelength. Note that the y-axis is in logarithm scale. The yellow and red curves denote the law for average temperature of the Sun and Earth. The blue line shows the Wien's displacement law explained below.
Wien's displacement law: explains the wavelengths over which the Plank's law is maximum as a function of temperature T:
The Stefan-Boltzmann law can be obtained by integrating the Planck's law in irradiance units () over all wavelengths or frequencies as follows:
As a result, one can see that the Stefan-Boltzmann law describes radiative energy as irradiance across all wavelengths, where is the Stefan-Boltzmann constant.
Graybody: Unlike a blackbody, natural objects are gray bodies, as they only absorb and thus emit a fraction of incident radiation and not all of it like a balck body. Radiation of a graybody is different than a blackbody up to a proportionality constant:
where the dimensionless coefficient is called spectral emissivity. As is clear, emissivity is an spectral quantity, which varies at different wavelengths (frequencies). Therefore, one can integrate spectral emissivity over all wavelengths to obtain the so-called broadband emissivity.
As a result, the Stefan-Boltzmann law for a graybody can be expressed as follows:
Kirchhoff's Law: For an object with constant temperature, the absorbed radiative energy is equal to the emitted radiative energy. Therefore, absorptivity and emissivity are equal , when an object is in thermodynamic equilibrium.

1-3 Short- and Long-wave Radiation

In study of the Earth systems, we typically divide the EM spectrum into shortwave and longwave segments. The main reason is that the temperature of the Sun and the Earth are drastically different and thus the wavelengths of their radiation. According to the Wien's displacement law we have:
Solar shortwaves () versus terrestrial longwaves ().
It can be shown that the solar radiation above 4 is almost negligible. Therefore, by convention, we have:

1-4 Absorption, Reflection and Transmission

Any incident radiation upon a medium will be absorbed, scattered, or transmitted.
The radiative energy balance at an interface can be written as:
where is absorptivity, is reflectivity and denotes transmissivity.
Often some of the terms in the above radiation balance equation are negligible:
Schematic of absorption and scattering in atmosphere and earth's surface

4- Solar Radition at Top of the Atmosphere

To understand hydrologic cycle, it is key to know that how much of the Sun's energy will reach to the earth's surface and what will be its distribution. Using the Stefan-Boltzmann law, we have:
In far-field radiation, basically it is reasonable to assume that the radiation beams are parallel. Therefore, given that:
Assuming that , one can easily compute the intensity of solar radiative energy at top of the atmosphere as follows:
Recent detailed calculation by NASA using satellite observations shows that the solar constant varies between 1361-1362 at the top of the atmosphere. However, it is still customary to use 1367 in hydrologic analysis as the solar constant.
Question: Why does the earth surface receive a non-uniform amount of radiative energy?
Schematic of global insolation and effects of earth geometry on distribution of the solar radiation on earth surface (left). The earth axis of rotation is tilted with respect to its orbital plane, which is in fact the main reason for seasonality of earth's climate (middle). The solar zenith angle (), encodes the reduction on the intensity of radiation over higher latitudes, while is the so-called solar azimuth angle, which determines the planar orientation of solar incident radiation with respect to north.
Insolation () is the intensity of Sun's radiative energy received at top of the atmosphere that accounts for the fact that the incident radiation beams are not always perpendicular to the surface:
Thus and the solar insolation is,
Summer Solstice (June 20 or 21) refers to the longest day in the Northern Hemisphere. As is evident, summer solstice occurs while the Earth Northern Hemisphere is at its closest distance to the Sun and receives its maximum annual daily radiation budget. Analogously, while the Northern Hemisphere is at its longest distance with the Sun throughout the year, it experiences the longest night of the year. The day of Dec 21 or 22 is the Winter Solstice in the Northern Hemisphere while it is summer solstice in the Southern Hemisphere.
When the Earth's celestial/equatorial plane passes through the center of the Sun, the day and night are equal. These days (Sep 22 or 23 or Mar 20 or 21) are called Equinoxes.

5- Net Radiation at the Earth Surface

In order to understand hydrologic fluxes at the surface, we need to know that how solar radiation partitions between the Earth's surface and its atmosphere. This partitioning is dynamic depending on the earth rotation and atmospheric and surface compositions such as cloud and snow cover. Recent analysis of the attenuation of the incoming solar radiative energy within the Earth's atmosphere shows that, on average, Earth has a planetary albedo of around 30%, while nearly 23% of solar insolation is absorbed by the atmosphere and clouds, leaving around 47% of the insolation reaching the Earth's surface (Trenberth, 2009).
Estimates of planetary radiation budget (Trenberth, 2009).

5-1- Net Radiation

More specifically, we must determine the net input/output of short and longwave radiation at the Earth's surface, known as the Net Radiation .
where is the net radiation at he Earth's surface, is the incoming shortwave radiation, is the earth surface albedo, is the longwave surface emissivity, is the incoming longwave radiation, and is the outgoing longwave radiation.
There are many empirical relations to estimate based on location, season, and atmospheric condition, while the longwave components can be estimated using the Stefan-Boltzmann law
where and represent the average air and surface temperatures in Kelvin and is the Stefan-Boltzmann constant.
The longwave surface emissivity for different land surface type is as follows (Aria 2001):
The surface albedos for different land surface type is as follows (Aria 2001):
Longwave atmospheric emissivity depends on profile of atmospheric consituents and temperature. However, the most important variable is water vapor concentration. There are multiple empirical models for . Some of them are noted below:
Example Problem 3.1: For a vapor pressure of 1250 [Pa], temperature of 25 [C] and air pressure of 98000 [Pa] estimate the clear-sky downwelling longwave radiation using the Brutsaert and Brunt emissivity models. How do the estimates differ and how would they change in cloudy conditions? Also compute the upwelling longwave radiation from the surface assuming it is covered by a long-grass surface with a surface temperature of 28 degree Celsius.
Solution:
% Brutsaert model
e_a = 1250/100; % 12.5 [mbar] of air vapor pressure
T_a = 25+273.15; % [K]
eps_0 = 0.90; % surface missivity [-]
T_g = 28+273.15; % ground temperature [K]
sig = 5.67*10^-8; % Stefan-Boltzmann constant
eps_a_brut = 1.24*(e_a/T_a)^0.14;
Rl_brut = eps_a_brut*sig*(T_a)^4; % [Wm^-2]
Rl_brut
Rl_brut = 356.3613
% Brunt model
eps_a_brun = 0.605+0.048*(e_a)^0.5;
Rl_brun = eps_a_brun*sig*(T_a)^4; % [Wm^-2]
Rl_brun
Rl_brun = 347.1035
% Long-wave radiation increases for a cloudy sky as the concentration of
% water vapor and liquid water increase.
% Upweling longwave surface radiation:
R_ul = eps_0*sig*T_g^4; % [Wm^-2]
R_ul
R_ul = 419.7175

5-2-2 Surface Energy Balance

The net radiation at the earth surface is partitioned into three components to create the Surface Energy Balance (SEB):
where is the latent heat flux, H is the sensible heat flux, and G is the ground heat flux. We will discuss these very important components in the lessons on evaporation. However, as way of review, below is a figure of the estimated global energy budget containing the components discussed in this section.
The global annual mean Earth's energy budget for 2000 to 2005 (Trenberth et al. 2009).
Example Problem 3.2: We have the following informationL vapor pressure of 1500 hPa, air temperature is equal to , P = 98,000 Pa. The surface soil tempertaure is . Estimate the clear sky downwelling long-wave radiation assuming that , where is the partial water vapor pressure in mbar and air tempertaure is in Kelvin. Additionally, compute the outgoing long-wave radiation with earth surface emissivity . Compute the net radiation when the shortwave incoming solar radiatoin is 410 and the surface albedo is 0.15. Is the surface soil warming or cooling?
Solution:
clear
R_is = 410; % [W.m^-2]
e_a = 15; % [mbar]
T_a = 20; % [C]
T_0 = 30; % [C]
sig = 5.67*10^-8; % [W.m−2.K−4]
eps_0 = 0.98; % [-]
alpha = 0.15; % [-]
eps_a = 1.24*(e_a/(T_a+273.15))^0.14; % air emissivity [-]
R_il = eps_a*sig*(T_a+273.15)^4; % incoming long-wave radiation
R_ol = eps_a*sig*(T_0+273.15)^4; % outgoing long-wave radiation
R_n = (1-alpha)*R_is+eps_0*R_il-R_ol % Net radiaton at the surface